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Abstract 

We show that an apphed charge current in a magnetic nanowire containing domain wahs (DWs) 
results in an interaction between DWs mediated by spin-dependent interferences of the scattered 
carriers. The energy and torque associated with this interaction show an oscihatory behaviour as 
a function of the mutual DWs orientations and separations, thus affecting the DWs' arrangements 
and shapes. Based on the derived DWs interaction energy and torque we calculate DW dynamics 
and uncover potential applications of interacting DWs as a tunable nano-mechanical oscillator. We 
also discuss the effect of impurities on the DW interaction. 

PACS numbers: 75.60.Ch,75.75.+a,75.60.Jk 
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I. INTRODUCTION 



Domain walls (DWs), i.e. regions of noncollinearity separating areas of different homoge- 
nous magnetization directions, are important from a fundamental and application point of 
viewii^i^. This is particularly the case at low dimensions, as in magnetic nanowires carriers 
turn out to couple strongly with DWs^ leading to a marked influence on the wire's transport 
properties, e.g. DW magnetoresistances in the range of 1000% were reported^'^'^ . As this 
coupling is associated with a change of the carriers' spin, it results in a current-induced spin 
torque acting on the DW and consequently in a current- induced DW motionii^. Based on 
these facts magnetic nanowires with a series of DWs can be utilized as a "racetrack DW 
memory"-. The DWs' motion is current-controlled; DWs separated by rather small distances 
are addressable thus allowing for a high memory density. 

In another context it is established that strong carrier scattering and interference results 
in long-range interactions between impurities on metal surfaces. This interaction governs 
the impurities geometric arrangements and growth^'^i^'^i'^'^. The question of whether 
and how the carriers' spin dependent scattering mediates interactions between DWs is still 
outstanding and should be addressed here. Clearly, the answer is of vital importance for 
a high-density nanowire-based racetrack memory and adds a new twist on interference- 
mediated interactions. We focus on the current-induced part of the coupling between 
neighboring DWs in a magnetic nanowire. Based on our results we identify the following 
mechanism of the DWs coupling: Upon scattering from the first DW a carrier spiral spin 
density builds up. This acts as a spatially non-uniform torque on the second DW whose 
energetically stable shape and position show therefore a non-uniform dependence on the 
distance from the first DW. This is different from the spin-torque transfer in bulk spin 
valve systemsi^ii^ii^ or magnetic tunnel junctionsi^ii^ insofar as in our case the DWs spatial 
arrangement, in addition to the magnetization direction, is current controlled^. We develop 
a theoretical framework to calculate the DWs current- induced effective potential and find it 
oscillates with the distance of DWs and their mutual polarization directions. This interaction 
we employ to study the DWs dynamics. As an application we propose the use of this new 
effect as a tunable, current-driven two-DW magnetic nano-oscillator— i^SiiSi; with a radiation 
emission dependent on the DWs positions in the various possible stable configurations. 
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II. THEORETICAL MODEL 



We consider a magnetic nanowire with two DWs when an electric current / is transmitted 
through the wire (a schematic is shown in Fig. [T]). When the distance Zq between DWs is 
larger than the phase coherence length L^, DWs act as independent scatterers. For zq < 
the current transmission mediates DW coupling. For definiteness, we assume that one of 
the DWs (located at z = 0) is pinned, e.g. by a geometric constriction, and concentrate 
on the effect of the current on the second DW, initially (i.e., for 1 = 0) located a.t z = zq. 
For / = each DW has an extension L. The transverse dimensions of the wire should 
be smaller than the exchange length and the Fermi wave length of the carriers, a situation 
realizable for magnetic semiconductors. The Hamiltonian H of independent carriers coupled 
(with a coupling constant J) to a spatially non-uniform magnetization (DWs) profile NL[z) 
is modeled by (we use units with h = 1) 



H = J dz al^{z) 



Saf^- J ■ M(2) 

2m 



ap{z), (1) 



where aj^ and aa are the creation and annihilation operators of electrons with spin a. Ap- 
plying a local gauge transformation T(2)^>2^>2i -we obtain instead of the nonuniform magne- 
tization a Zeeman splitting term and a spin-dependent spatially varying potential Ua/3{z), 
which for kpL > 1 can be treated perturbatively^i^^i^iiS^ (^kp is the Fermi wave vector). 
Note, for a sharp domain wall, i.e. for kpL < 1, the formalism of Ref.— can be adopted. 
T{z) is obtained from the requirement T'^iz) cr ■ n{z) T{z) = a^, where n is the unit vector 
along M, Nl{z) = Mn{z). The transformed Hamiltonian H = T'^lz) HT{z) reads 

H = j dzaliz) 

with the perturbation given by 



-^5^p + U^p{z)- JMalp ap{z) (2) 



U{z) = -^[2Ad, + {d^) + A^] , (3) 

and A{z) = T'^iz) VzT{z) is a gauge potential. For a wire with two DWs we parametrize 
the magnetization profile by the angles ^p{z) and 6{z) (cf. Fig{T]) 



11(2;) = ( cos 6^ sine/), sin 9 sin ip, cos ip), (4) 

T 



(p{z) = COS ^ ( tanh [— ] ) -|- COS ^(tanh[ — (5) 



(See reference^l and references therein.) The angle 6{z) describes the relative orientation 
between the wall pinned at 2; = and the other situated around z = zq. We set 61 to zero 
at the first wall and 62 = Oq around the second (see Fig. [T] ). For 60 = 11 the walls are 
antialigned. For zq < L DWs may merge, hence we consider the case zq > L for which we 
may write U{z) ~ Ui{z) + U2{z), where (j = 1,2) 
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H — smfe',-. (6) 

4m 2m J 

This approach is generalizable to any number of DWs, which are sufficiently far apart. As 
shown in Refs.— "^"^ for a single DW, for kpL > 1, i.e. when Nl{z) hardly varies within 
fcj^^(adiabatic DW), the terms in Eq. ([6]) proportional to fi{z) are negligibly small and 
a perturbative approach is appropriate for treating the electron scattering from the DWs 
potential (Eq. ([6]))^. Assuming ip^{z) to be the wave function of an independent electron 
with energy e in the wire without the DWs, we find the first-order correction due to the 
perturbation Ui{z), i.e. due to scattering from the first DWs, as 



dz' Ge{z,z')Ui{z')^^{z'). 



(7) 



The Green's function corresponds to the unperturbed Hamiltonian with U{z) = 0. It is 
diagonal in spin space with elements 



Geaiz, z') 



kcr 



(8) 



where ~ k'^ 



I m 
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for lifetimes To- ^ Ep^, and k^ = [2m (e + /i ± JM)]^^^. Hence 



dz' I "'-T ^^'^ 



(9) 



and 



dz' 



^^e'^^\^-''\^'Az)\e'^^^' 



i^ik^\z~z'\ [^/^(^)j2g%y 



(10) 



for incoming electrons of spin up and down, respectively. 

The interaction energy of the two DWs due to the single scattered state ip£a{z) = ip^^{z) + 
S'ip£fj{z) is calculated as 

POO 

AE^= I dz54,{z)U2{z)5iJea{z). (ll) 



Summing up the contributions of all scattering states in the energy range between ep and 
Ep + eA0/2, for an applied voltage eA0/2 ^ ^i?, we obtain the current-induced coupling of 
the DWs as 

V2it\ ft n 
where v^j = k^/m is the velocity of electrons at the Fermi level. 




III. NUMERICAL EXAMPLES 

Magnetic semiconductors^i^^ are most favorable for a sizable effect, for metallic wires the 
ID limit is also within reacb^. Here we use in the numerical calculations similar parameters 
as in RefQ, i.e. Xp = Q nm; a mean free path of / = 500 nm; an effective mass of m = O.Sm-e 
{rrie is free electron mass); L = Xp', JM = 15 meV; ep = 83.7 meV; and eA0 = O.lei;'.— 
The width of the wall may well be on the atomic size in the presence of constrictions^i^i^, 
i.e. well below the DW lengths in bulk materials. In such a situation, the DW-interaction 
increases due to the strongly enhanced DW scattering2^i^i2ii2^. The interaction energy 
Figure [1] depends periodically on the DWs mutual angle 6*0 and distance Zq, which results in 
an oscillating motion of the DW along the axis z as well as an oscillating direction of DW 
polarization. 

Now we focus on the effect of DW scattering on the electron spin density, leading to a 
nonequilibrium spin accumulation and to a spin torque acting on the wall. Subsequently, 
we study the dynamics of the DW related to the DW coupling. 

The spin-density due to the single transmitted wave of spin a is 

S^{z)=^l{z)T{z)aT\z)ij,^{z), (13) 
and the total current-induced spin density is^ 



We find that the correction to the spin density follows the magnetization profile with addi- 
tional Friedel oscillations, which are a superposition of two waves with periods kp^ and /c^^J. 
The oscillations in the spin density are smaller in magnitude than the overall spin density 
profile and decay with increasing z. 
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FIG. 1: (Color online) Top panel: A schematics showing the DWs magnetization profile (thick 
arrows). L is the DW width, zq and 9o are respectively the DW position and orientation with 
respect to the DW at z = 0, / is the current direction. Lower panel: Interaction energy AE{zq, 6q) 
as a function of zq and ^o- Solid curve is for zq = 30 nm, the dashed is for zq = 37.5 nm, and the 
dotted is for zq = 45 nm. 

We calculate the current-induced torque acting on the second DW at z from 

AT(^, ^0, ^o) = -— M(z, zo, ^o) X AS(^, zq, ^o), (15) 

where 7 = gfiB, 9 is the Lande factor and /is is the Bohr magneton. We assumed a 
thin nanowire with a cross section of = 100 x 20 nm^ as in Ref.-. In Eq. (15) AS is 
the correction to the electron spin density due to scattering. The calculated torque on the 
second DW is shown in Figs. [2] and [3], where zo = 50 L and M ^ 5.56 x lO^Am"^ were used^^. 
The correction to the spin torque shows that the force upon the DW depends strongly on 
their relative polarizations. 

To inspect the current-induced dynamics of the DW at 2; = ^o, we evaluate the accumu- 
lated spin density that acts on the DW at z = zq. The DW magnetization dynamics are 
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FIG. 2: (Color online) The x-component of the current-induced spin torque, as defined in Eg. 1151 
acting at the second domain wall as a function of z and ^o- 
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FIG. 3: (Color online) The z-component of the current-induced spin torque around the second 
domain wall as a function of z and 9q. 

then modeled using the Landau-Lifshitz equation^ 

(9tM = -— M X S[M]. (16) 

O'cs 

As an initial condition we assume that the magnetization profile in the wire without electric 
current is described by Eq. ([5]). The results for the time dependence of the magnetization 
are shown in Fig. H] for the centre of the DW, z = z^. We should note that the relative 
orientation of the walls at the start of motion does play a role in the type of motion we see. 
Here we present it for an arbitrary configuration. As we move away from the centre of the 
DW, the relative orientation becomes increasingly irrelevant. At z = zq + L the motion is 
the same regardless of the value of ^o- 

Analyzing different magnetization components, we find their motions have different fre- 
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FIG. 4: (Color online) The time dependence of the magnetization with the initial condition for the 
second wall to be at an angle of = 7r/4 to the first wall. This is the solution to Eq. [T6l The 
solid curve is the x-component, dashed the y-component, and dotted the z-component. Taken at 
the centre of the domain wall. 

quencies and a different form. No lateral movement or permanent distortion of the DW is 
observed, we see only oscillations. At the edge of the DW wall there are small rotations of 
the magnetization, regardless of the initial conditions. If we move far from the DW then all 
the dynamics of the magnetization vanish. 

This is in contrast to the case where we do not include the first domain wall. In this case 
there is no motion near the centre of the domain wall at all. Furthermore the oscillations in 
the magnetization towards the edge of the domain wall are much slower (by several orders 
of magnitude) than exhibited here. 

Extending our analysis to include the effect of magnetic anisotropy we write 

dtM = -^M X S[M] + X xM,. (17) 

We take the anisotropy constant K' = —10, and therefore the x-axis as a hard magnetization 
axis. Figure [5] shows the effects of anisotropy on the domain wall motion. The anisotropy 
dampens motion in the x-direction, thus exacerbating the y and z oscillations. This is also 
in contrast to the case where we ignore the first domain wall. In this case, although the 
anisotropy does introduce motion around the centre of the domain wall it does not involve 
a decaying x-component, see Fig. [6l 
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FIG. 5: (Color online) Same as in Fig. U] when the anisotropy is included, see Eq. [T71 The solid 
curve is the x-component, dashed the y-component, and dotted the z-component. Taken at the 
centre of the domain wall. 
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FIG. 6: (Color online) The magnetization dynamics of the second DW without the influence of the 
first DW, when the anisotropy is included, see Eq. [T71 The solid curve is the x-component, dashed 
the y-component, and dotted the z-component. Taken at the centre of the domain wall. 

IV. SUMMARY 

A current through a magnetic nanowire containing DWs results in a DW interaction medi- 
ated by the scattered charge carriers. We developed a method for calculating the interaction 
energy and the consequences of this new coupling mechanism. The DWs interaction energy 
oscillates as a function of the DWs mutual orientation and distance. This has immediate 
consequences on how DWs rearrange upon applying a bias voltage and on the fundamental 
limit of the DWs packing density. In fact, different parts of the DW oscillate at different 
rates and in different ways: becoming more regular, smaller, and quicker away from the DW 
centre. The nonequilibrium DWs oscillations around the energy minima generates radiation 
with a frequency dependent on the applied bias voltage, DW length and scattering strength. 
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These parameters are externally tunable for utilizing the interacting DWs as a versatile 
radiation source. For an experimental realization magnetic semiconductors^i2^i21i2^i2£dS are 
favorable, our results are in the range already achievable^. Extension to the metallic case is 
straightforward, though DW lengths may not be easily fabricated on the required scale. In 
this case the results remain qualitatively similar. Furthermore anisotropy will completely 
dampen any DW oscillations as motion in a plane becomes much harder due to the much 
larger magnetization size, M = 1.72 x lO^Am"^ for Fe. 

In our numerical simulations we used parameters of a magnetic semiconductor with rela- 
tively large electron wavelength, X < L, and much longer mean free path / ^ A. The latter 
can be realized in case of small density of impurities and defects. However, magnetic semi- 
conductors like GaMnAs are usually strongly disordered, and instead of a strong inequality 
one may find / > A. In this case, the phase of the current-induced spin density wave at 
a distance zq > I will be affected by impurities, and, therefore, one can expect that the 
disorder- averaged interaction between two DWs at a distance zo is suppressed by the factor 
g-2o/«_ However, we should stress that the real interaction between two DWs depends only 
on a given realization of the disorder and therefore is not damped by the impurities. This 
effect is analogous to the nondamping of the RKKY interaction between magnetic impurities 
in disordered metals^i^. The detailed analysis of this phenomenon is beyond the scope of 
this paper. 
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